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The problem of antiproton scattering on the molecular hydrogen is investigated by means of wave packet 
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work. Excitation and dissociation probabilities of the molecular hydrogen for collision energies in the ultra 
low energy regime below 10 eV are computed. 
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I. INTRODUCTION 

The antimatter-matter asymmetry is an active area of 
physics adressing fundamental questions concerning the 
evolution of the universe, the validity of the CPT sym¬ 
metry and gravitational differences between matter and 
antimatter^. 

Even though our universe seems to be solely composed 
of regular matter it is possible to create and study anti¬ 
matter in the laborator}!^. The recent advances in the 
field of antimatter include the production and trapping 
of antihydrogeiP as well as cooling to a stable electronic 
stat^. Antihydrogen is one of the simplest examples 
of an atomic system composed of antimatter. However, 
to increase the understanding of antimatter systems the 
interaction between matter and antimatter is of special 
interest. Mixed matter-antimatter many-body systems 
provide a field of study, which exceeds the pure annihila¬ 
tion effects and allows for comparison with the dynamics 
of ordinary molecular systems. Examples of such many 
body systems includes, e.g., antihydrogen- hydrogerP, an¬ 
tiprotonic heliunJ^ and various antiproton - gas colli¬ 
sions investigated by the ASACUSA collaboration in the 
antiproton collisional kinetic energy regime of 2-11 keV^. 

In this work we have investigated the collision an an¬ 
tiproton (p) with molecular hydrogen as a model sys¬ 
tem for more complex antiparticle - molecular interac¬ 
tions. The H 2 -P system is interesting from a practical 
point of view as residual H 2 gas exists in most high en¬ 
ergy colliders, and antiprotons are created in various ex¬ 
periment^ resulting in scattering of antiprotons on H 2 
molecule^. Previous work on this system focused on the 
high-energ}ffll and the low-energy collision region investi¬ 
gating antiproton and proton cross sections for ionization 
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and excitation of hydrogen m olecul es as well as energy 
spectra of the ionized electron d^^ l ^^f ^. 

In previous work the influence of collisions in the en¬ 
ergy regime above 100 eV was investigated which mainly 
leads to ionized reaction product^. However, the ultra 
low energy regime below 10 eV^ opens up a new regime 
of matter-antimatter chemistry, since this energy is com¬ 
parable to the chemical bond energy. In this region the 
molecular hydrogen is vibrationally excited, but remains 
in its electronic ground state. The antiproton behaves 
like another nucleus in the molecular system, whereas its 
negative charge produces an interaction potential, which 
is on the order of magnitude of an electronic state in an 
atom. 

In this paper we have computed potential energy sur¬ 
faces (PES) for the lowest electronic states of the H 2 -P 
system and performed wave packet simulations of the col¬ 
lision on the electronic ground state. The wave packet 
simulation provide information about probability distri¬ 
butions of the vibrational states of the H 2 molecule after 
the collision, as well as the dissociation probability. The 
paper is organized as follows: In section |TI| the quan¬ 
tum dynamical model is introduced. In section HI the 
potential energy surfaces (PES) of the H 2 -P system and 
the results of the wave packet calculations are presented 


followed by a discussion in section IV 


II. QUANTUM DYNAMICS MODEL 

We have considered a T-shaped configuration of the 
H 2 -P system corresponding to a zero impact parameter, 
as illustrated in Eig. which is expected to have a max¬ 
imum interaction between the antiproton and the vibra¬ 
tional degree of freedom of the molecular hydrogen. Ja¬ 
cobi coordinates are used to describe the inter-atomic 
distances between the H-atoms {Rhh) and the antipro¬ 
ton distance to the center of mass of H 2 {Rp)- The angle 
of the H-H bond axis and Rp is fixed to be 90 degrees. 
This choice of coordinates focuses specifically on the vi¬ 
brational and dissociation dynamics of H 2 in collision 
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FIG. 1. Jacobi coordinates for the H 2 — p system in a T- 
shaped configuration. The internuclear separation between 
the hydrogen atoms Hi and H 2 are varied along the Rhh 
axis and the distance from the barycenter of the two hydrogen 
atoms to the antiproton is the Rp axis. 


event with p. Other reaction pathways like the proto- 
nium format iorP^ are neglected in this model. 

To simulate the collision between the hydrogen 
molecule and the antiproton we solve the time-dependent 
Schrodinger equation for the nuclei 

Rp, t) = H^{Rhh, Rp, t ), (1) 

numerically. Atomic units are used (i.e. h = rrie = e = 
47reo = 1) in what follows. The Hamiltonian in Jacobi 
coordinates reads: 

H = Trp + Trhh + ^ {Rhh, Rp) = (2) 

15^ 1 ^ 

“ ~ + ’ 

where Rhh, and Rp are the Jacobi coordinates as defined 
in Fig. The reduced masses and pp and phh are given 
by 


Trip 

I^hh = ^ 

2ml 2 

^ - 7 = 77^^ 

^ {rup + 2mp) 3 


(3) 

(4) 


where nip is the rest mass of the proton. The choice 
of the Jacobi coordinates ensures that the system is in 
the center of mass frame. The molecular potential en¬ 
ergy surface (PES) V represents the adiabatic electronic 
ground state potential of the system and is obtained by 
solving the electronic Schrodinger equation under the as¬ 
sumption that the Born-Oppenheimer approximation is 
valid. The collision of the molecular hydrogen with the 
antiproton in the given coordinate system Rhh, Rp is 
simulated by solving eq. numerically with a Cheby- 
shev time propagation schem^I^^ A perfectly matched 
layer (PML)P^ at the edges of the numerical grid has 
been used absorb the wavefunction and to avoid spuri¬ 
ous reflections. 


The initial condition is chosen to meet the given phys¬ 
ical conditions of an antiproton approaching a hydrogen 
molecule. The initial wave function can be written as 
a product state of the H 2 vibrational ground state wave 
function and a Gaussian wave packet representing the 
antiproton: 


^{Rhh,Rp, t = 0) = (5) 

y (jy2ii 

The collision energy po = 1‘^Rp enters eq. [^through 

the momentum p = hk. The width a of the wave packet 
is assumed to be 0.3 A, which corresponds to a width of 
35meV in kinetic energy. In the limit of i?o ^ 00 i.e. 
of a separated reactants the PES becomes flat and the 
Schrodinger equation becomes separable in the molecular 
coordinates justifying the choice of the wave-packet, eq. 
1^ A value of i?o = 10 A is chosen as a good approxima¬ 
tion on a finite surface used in the numerical model. 

By describing the system with the non-relavistic 
Schrodinger equation, the particle anti-particle annihila¬ 
tion of the proton and the antiproton is neglected. How¬ 
ever, a direct p-p collision has a low probability, specially 
for the given setting (impact parameter of zero). This is 
justified by the wave packet simulations which show that 
the probability density of the wave function is low at 
the positions representing a direct p-p contact. A com¬ 
parison of the simulation time scale for a single collision 
event with experimentally determined annihilation time 
constant^l^ shows that the annihilation probability is in¬ 
deed low and annihilation effects are expected to be small 
and to have no significant influence here. 


III. RESULTS 
A. Potential energy surfaces 

The electronic ground state and the first five excited 
states potential energy surfaces have been calculated with 
the full configuration interaction (ECI) method using 
MOLCAS program packag^^. The Atomic Natural Or¬ 
bital basis set of valence triple zeta accuracy with po¬ 
larization functions (ANO-L-VTZP) basis sefP^ is used 
to describe the two electrons. The influence of the an¬ 
tiproton on the electronic states of the H 2 molecule is 
simulated by adding a negative point charge to the geom¬ 
etry definition. A comparison between the ground state 
and the first few excited electronic state potential energy 
surfaces (PES) is presented in one-dimensional slices at 
different antiproton distances in Eig. i In Fig. [g one can 
see that the electronic ground state is well separated from 
the excited state surfaces and no avoided crossings or con¬ 
ical intersections between the ground state and an excited 
state occurs. However, curve crossings can be observed 
between the 1^11^ and the 2^n^ state at Rhh ~ 1.6 A 
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FIG. 2. Slices through the PES for different values of Rp. (a) 0.8 A, (b) 2.0 A, and (c) 6.0 A. The ground state (l^Eg) and the 
first - to fifth excited states (H 2 state symmetries l^Eu, 2^X]g, 2^11^ and 2^Eg) energy are varied along the internuclear 

separation between the hydrogen atoms Rhh- The labels of the electronic states refer to the electronic states of unperturbed 
H 2 in the Dooh point group. A potential energy of OeV refers to the limit of separated particles. 


and 2^n^ and the 2^E^ state at Rhh ~ 1.2 A. Note that 
for the hydrogen molecule only avoided crossings can oc¬ 
cur, wherea s the triatomic H 2 -P system can have conical 
inter sect ionPSmi between states of the same symmetry. 
As the antiproton approaches the H 2 molecule the PES 
energy and the difference between the vertical excitation 
energy and the groundstate energy decreases, while the 
l^n^, 2^n^ and 2^ Eg states increase their separation in 
energy. For the 2^ Eg and excited state crossings 

have been observed. The groundstate energy curve be¬ 
comes identical to the ground state curve for unperturbed 
H 2 as ^ 00 . The dissociative limit of the H 2 ground- 
state approaches -27.21 eV, which corresponds to the sum 
of two H atoms in their groundstate, as Rhh c>c. 

A two-dimensional plot of the electronic ground state 
PES in the T-shaped configuration, is shown in Fig. A 
deep potential well is found where the antiproton passes 
through the axis of the molecular hydrogen caused by 
an attractive Coulomb interaction between the protons 
and the antiproton. The theoretical value in the origin of 
the coordinate system {Rhh^ Rp} can be directly derived 
from the Coulomb interactions Vc by taking the limit 
{Rrh^Rp} 0+ for the interaction potential: 



FIG. 3. PES for the H 2 —p systems ground state, calculated in 
dependence of the H 2 separation and the antiproton distance. 
For clarity the potential is truncated above +50 eV. Note that 
the potential converges to —00 at (0,0). 




lim Vc = lim 

{Rhh, Rp}^0+ {rip, r2p, ri2}^0+ Tip 


+ 


CpCp ^ G-p 

r2p ri2 


1 and the antiproton, proton 2 and the antiproton, and 
^ the distance between proton 1 and 2 respectively. Here 

_ _ ^ Cp = — 1 are the charge of the proton and the 

{Rhh,Rp}^o+ Rhh \/R- + {Rhh/‘^)‘^ antiproton respectively. From eq. [^it becomes clear that 
^ the limit diverges to minus infinity. This is indicated in 

where rip, r 2 p and ri 2 are the distances between proton Fig. [^showing a potential hole diverging to minus — 00 . 
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B. Quantum Dynamics Simulations 

The H 2 antiproton collision is simulated by solving 
eq. □ numerically using the Chebyshev time propaga¬ 
tion schem^i^, combined with the Fourier methocP^l, for 
different collision energies. The initial conditions are con¬ 
structed from eq. [^and cover collisions with an incoming 
energy between 0.01 and 10 eV. The PES which enters 
through the Hamiltonian (eq. is the l^Ei electronic 
ground state as it is shown in Fig. We assume that 
the electronic ground state of the H 2 -P system is sufficient 
to describe the collision in the energy regime considered. 
This is justified by the absence of curve crossings between 
the ground state and the excited states. 

The time evolution of a wave packet with an initial 
kinetic energy of 3eV is shown in Fig. [^for 3 different 
times. As the wave packet approaches the deep potential 
well the strong attraction of the Coulomb force becomes 
visible (Fig. &))• When the antiproton transverses the 
potential singularity, the wave function is either reflected, 
transmitted, or deflected towards the H 2 dissociation. As 
the potential is a conservative potential in this model the 
wave packet can not be trapped in the location around 
the potential singularity. Eventually it will be scattered 
into separate fragments. The quantitative results of the 
inelastic scattering event will be presented in the follow¬ 
ing sections. 


1. H 2 vibrational excitation 

The vibrational excitation after the collision event is 
analyzed separately for the reflected and the transmitted 
part of the wave function. In the asymptotic limit of a 
large |i?p| the Schrodinger equation for the Hamiltonian 
(eq. becomes separable and the vibrational excitation 
can be expressed in terms of vibrational eigenfunctions 
(I>u{Rhh) of H 2 . The probability of finding the system 
in a vibrational eigenstate at a certain time t is given 
by the projection against the total nuclear wave function 

= J\(l^i^{RHH)'^{RHH,RpoR)fdRHH- (7) 

The distance from the H 2 center is chosen to be Rpo = 10 
A, which is far enough from the interaction region such 
that the separability assumption is sufficiently fulfilled. 
The total probability is then obtained by integrating 
eq. [^over tim^^ 

= \c,,{t)\‘^dt. (8) 

The time end point r (500 fs) is chosen such that less 
than 1% of the wave packet remains after absorption at 
the PML. A time step 242 as has been chosen and is 
sufficiently short to allow for an accurate numerical inte¬ 
gration over time. 


The antiproton is scattered in all directions and inte¬ 
gral bounds have to set such that the reflected part of the 
wave function and the transmitted part of the wave 
function Tt is captured. The sum over the coefficients 
corresponding to the these probability distributions is the 
normalization coefficient we use, i.e. 

Ar^ = y](/ \c^{t)\^dt+ f \cl{t)\^dt) . (9) 

\Jt=o Jt=o J 


The normalization is chosen such that the sum of 
the probabilities for the reflected and transmitted part is 
unity. 

The total probability for the elastic scattering P^{h' = 
0) is given by 




which, together with the dissociation probability Pdiss 
described in section IIIIB2I gives the resulting inelastic 
scattering probability. 

The vibrational distributions for the transmitted and 
reflected part are shown in Fig. Until 0.1 eV elastic 
scattering dominates the process and H 2 remains in its vi¬ 
brational ground state. At 0.28 eV the reflected part be¬ 
gins to slightly broaden its vibrational distribution while 
the transmitted part still remains in z/ = 0. This is the 
main difference between the reflected and the transmit¬ 
ted vibrational distributions, which otherwise are very 
similar. As the kinetic energy is increased further, u dis¬ 
plays a broad distribution of the vibrational states. For 
comparison the mean vibrational quantum numbers 


(z/) = vP^ 


u=0 


( 11 ) 


of the H 2 vibrations after the collision are shown in Fig. 
m Due to the small difference between the reflected and 
transmitted parts only the sum of both is shown. The 
mean vibrational quantum number increases linearly un¬ 
til the H 2 dissociation energy of 4.52 eV and peaks around 
7 eV with a maximum around (u) = 10. The mean vi¬ 
brational quantum number can be regarded as a good 
measure until about 6eV where the dissociation proba¬ 
bility reaches a threshold and increases. 


2. H 2 Dissociation 


The dissociation probability of H 2 is analyzed by 
using the probability flux derived from the continuity 
equation^ 


dt 


+ V ■ i = 0 , 


( 12 ) 
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FIG. 4. Time evolution of the antiproton wave packet with a kinetic energy of po = 3eV at times t = 60 fs, t = 110 fs and 
t = 160 fs (from left to right). 



FIG. 6. The mean vibrational quantum number (z/) distribu- 
(a) tion versus initial antiproton kinetic energy po for the sum of 

the reflected and the transmitted wavefunction. 



(b) 


FIG. 5. The probability distribution for the reflected (a) and 
the transmitted (b) wavefunction for the H 2 molecule to end 
up in vibrational excitation state to 14 versus initial an¬ 
tiproton kinetic energy po. The probabilities of both diagrams 
sum up to one for a specific energy value. 


where j is the probability flux. Eventually, the probabil¬ 
ity flux F(t) reads: 

= J ^ {^{RHH,Rp,t)*f<i>{RHH,Rp,t)) dV, 

(13) 

where ^ denotes the imaginary part, and V is the vol¬ 
ume corresponding to a bound H 2 in the coordinate sys¬ 
tem Rhhi The operator T is the kinetic operator 
from the Hamiltonian (eq. |^. By using a step function 
Q{Rhh — Rhh,o) which drops to zero at Rhh = 15 A, 
the probability flux describing the dissociation of H 2 is 
obtained: 


2 

Rdiss{t) = - V / &{Rhh — Rhh,o) (14) 

^ J — 00 

X ^ (^^{Rhh, RpRTT^iRHH, RpR)^ dK. 

The integration over time yields the dissociation proba¬ 
bility: 


diss 



Rdiss {t)dt 


(15) 
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FIG. 7. Dissociation probability P(D) of the H 2 molecule ver¬ 
sus initial antiproton kinetic energy po, displaying a rapidly 
growing behavior after 6 eV. 


The dissociation probability is shown in Fig. It is close 
to zero until it reaches a dissociation threshold at 6eV. 
Between 7 and 10 eV the dissociation probability increase 
linearly to 25 %. Note that the dissociation energy of H 2 
is 4.52 eV. 


IV. DISCUSSION 

The ground state of the H 2 -P system is separated by 
more than 5eV from the first excited state, which can 
be seen from Fig. and thus the effect of non-adiabatic 
couplings can be neglected if the H 2 is initially in its 
electronic ground state. We notice, however, that the ex¬ 
cited states 2^n^ and 2^11^^ of the coupled system show 
clearly a region of curve crossing (Fig. [^b)), which needs 
to be taken in consideration for higher collision ener¬ 
gies. In a three body system an avoided crossing can 
become a conical intersection (CI]P^with a truly degen¬ 
erate point. The third degree of freedom necessary to 
form a Cl originates from an anti matter particle with 
a negative charge instead of an ordinary nucleus, mak¬ 
ing it a special case. This can also be viewed as a more 
general case of CIs. Other examples for introducing CIs 
through a non-standard nuclear degree of freedom are 
light induced CI^. 

The ground state potential in two dimensions (Fig. 
demonstrates the unique features of the interaction 
with a negative charge and the hydrogen molecule. The 
PES is characterized by its deep Coulomb hole diverg¬ 
ing to minus infinity. The form of the potential suggests 
that a bound vibrational ground state could exist. How¬ 
ever, this ground state - representing a captured antipro¬ 
ton - could be subject to strong annihilation effects and 
short lived. In contrast to a tri-atomic molecule the H 2 -P 
ground state is expected to be several keV deepP^ due 
to the mass of the antiproton (compared to an electron). 
Its determination and investigation however is beyond 
the scope of this paper. 


The analysis of the vibrational distribution of H 2 (Fig. 
shows a clear threshold in collision energy below 
which the reaction products are mainly in the vibrational 
ground state. This threshold value is between 0.28 and 
0.63 eV and is consistent with the the first vibrational ex¬ 
citation of H 2 at 0.55 e\^. However, the distribution has 
a slight asymmetry between the reflected and transmit¬ 
ted part at 0.28 eV. This indicates that the slightly longer 
interaction time for the reflected part allows coupling of 
some of the antiproton’s kinetic energy into the H 2 mode. 
For collision energies beyond 3.4 eV the vibrational distri¬ 
bution becomes more spread out, which can be explained 
by the anharmonicity of the H 2 : at higher collision en¬ 
ergy more vibrational levels become accessible. Here the 
mean vibrational quantum number (u) (Fig. provides 
a more clear picture. The increase in (z/) is close to a 
linear behavior for collision energies below 4eV. A com¬ 
parison with the energy of u = 7 (2.9 eV) and the col¬ 
lision energy leads to the conclusion that approximately 
90% of the antiprotons kinetic energy is converted into 
vibrational energy. 

At a collision energy of ~ 6eV the H 2 begins to dis¬ 
sociate (Fig. [^. However at least 8eV is necessary to 
dissociate about 10 % of the H 2 . Compared with the dis¬ 
sociation energy of 4.52 eV of H 2 this is a remarkable 
excess of energy which is needed to induce the H 2 disso¬ 
ciation. This indicates that the coupling between the an¬ 
tiproton mode and the H 2 mode becomes more inefficient 
with increasing collision energy. Such effects have been 
observed in other reactive scattering systems in the form 
of a dynamical barrieil^. A simple explanation for this 
phenomena could be provided by the fact, that the inter¬ 
action time between the H 2 and the antiproton becomes 
shorter at higher velocities and reduces the effective time 
to transfer energy between the modes. 


V. CONCLUSIONS 

To describe the dynamics of an antiproton colliding 
with a hydrogen molecule the electronic PES of the 
ground state as well as for the five lowest excited states 
were first computed. The quantum chemical calculations 
show that it is enough to include the ground state PES 
to describe the dynamics in the ultra low energy regime 
considered here. For higher collision energies it might 
be necessary to include non-adiabatic couplings, as indi¬ 
cated by the curve crossings between the excited states. 
The possibilities of finding conical intersection instead 
of avoided crossings, compared to bare H 2 , is an excit¬ 
ing example of non-conventional conical intersections. A 
thorough analysis will be conducted in a future paper. 
The results of the wave packet dynamics give insight into 
the energy redistribution during the collision. For ener¬ 
gies below the dissociation threshold the scattering event 
leaves a clear signature in the vibrational distribution of 
the H 2 system. The vibrational excitations are in accor¬ 
dance with the H 2 vibrational states. The system has a 
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clear threshold with respect to dissociation. The theo¬ 
retical calculations show that this threshold is 1.5 eV 
above the dissociation limit of H 2 . 
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